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Abstract 



■ We study the Nordstrom- Vlasov system which describes the dynamics of a self-gravitating 

ensemble of collisionless particles in the framework of the Nordstrom scalar theory of gravi- 
tation. If the speed of light c is considered as a parameter, it is known that in the Newtonian 
limit c — > oo the Vlasov-Poisson system is obtained. In this paper we determine a higher ap- 
proximation and establish a pointwise error estimate of order 0(c -4 ). Such an approximation 
q is usually called a 1.5 post-Newtonian approximation. 



1 Introduction and Main Results 



In astrophysics stars of a galaxy are often modeled by a collisionless gas interacting only by the 
gravitational fields which they create collectively; together this leads to the Vlasov-Einstein system. 
Recently, interest has arisen in a simplified but still relativistic model essentially going back to 
Nordstrom, see [T^j, in which the dynamic of the matter is coupled to a scalar theory of gravitation, 
the metric tensor potential is replaced by a scalar function and Einstein's equations are replaced 
by a wave equation. In 4J a reformulated version of Nordstrom's theory is presented in which the 
system reads 

S(f) - \s(<f>)p + 1C 2 V X <f>] • V p / = 4S(cj>)f ( 

(VNc) 




In the previous equations, / = f(t,x,p) gives the probability density to find a particle at time t 
at position x with momentum p where t 6 R, 16 M. 3 , p G M 3 . <p = <p(t, x) is the mean Nordstrom 
gravitational potential generated by the particles, c the speed of light. Moreover 

P 2 = \p\ 2 , 7 = (1 + c 2 p 2 )~ 1/2 , p = 7p, and S = d t + p ■ V x 
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is the relativistic free-streaming operator and p the relativistic velocity associated to the momentum 
p. 

Together with the initial conditions 



(jVNc)) is the Cauchy problem of the Vlasov-Nordstrom system. For a physical interpretation and 
a derivation of this system see [3]. 

In this formulation (jVNcJl exhibits many similarities to the relativistic Vlasov-Maxwell system, 
which also consists of a, in that case homogenous, Vlasov equation coupled to a linear hyperbolic 
equation. Thus, many techniques developed for the Vlasov-Maxwell system also apply for the 
Vlasov-Nordstrom system providing existence and uniqueness of local classical solutions as well as 
a continuation criterion, which guarantees that the solution is global in time if a control for the 
velocity support of the matter density is established, see [3] . Furthermore, there are results about 
existence of global weak solutions, see |H] , existence of global classical solutions in the 2D case, see 
|lf)j . for spherically symmetric initial data, see pQ, and for small initial data, see |0j. Concerning 
the Vlasov-Maxwell system we refer the reader to [7] and the references cited therein. On the 
other hand the same important pieces of the basic theory are missing, namely global existence of 
classical solutions with unrestricted initial data and uniqueness of weak solutions. However, it is 
remarkable that there is a blow-up result for ()VJNc|) if the sign in the force-term of the Vlasov 
equation is changed, see jl]. 

This paper is concerned with the non-relativistic limit of (|VNc|) . i.e. the limit c — > oo: From 
physical intuition the Vlasov-Poisson system is the non-relativistic limit of the Vlasov-Einstein 
system, a statement which has been made rigorous under certain circumstances in j22J. In [3| 
it has been shown that as c — > oo also solutions of (jVNcj) converge to a solution of a Vlasov- 
Poisson equation in a pointwise sense obeying an error estimate of order 0(c~ l ). Thus, (jVNc|) can 
be regarded as another relativistic generalization of the Vlasov-Poisson system. The situation in 
the Maxwell case is similar as solutions of the Vlasov-Maxwell system approach a solution of the 
Vlasov-Poisson system (in the plasma case, which differs from the Vlasov-Poisson system for the 
gravitational case by the sign in the force-term) in the same sense, see |23j . 

It is the goal of this paper to replace the Vlasov-Poisson system, the classical or Newtonian limit, 
by another effective equation to achieve higher order convergence and a more precise approximation 
(in the plasma case this has been done in [2]). This will lead to an effective system whose solution 
stay as close as C(c -4 ) to a solution of (|VNc|) if the initial data are matched appropriately. Thus, 
it is the 1.5 Post-Newtonian approximation, note that in the context of General Relativity the 
approximations are usually counted in powers of c~ 2 . 

In order to derive effective equations in the limit of small initial velocities we shall formally 
expand all quantities in powers of cr x 



f(0,x,p) = f°(x,p), 0(0, x) = 0°(x), d t cf>(0,x) = <j>\x) 



(1.1) 



/ = fo + c"Vi + c~ 2 h + ■ ■ ■ 

/I — /J, + c~Vi + C _2 /i2 H 

= 00 + C~Vl + C~ 2 2 H 



put this ansatz into (|VNc|) and derive an equation in every order of c 
Therefore, we obtain in the orders c 2 and c 1 



,-i 



-A X O = O, - 



A, 0x = 0, 
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thus, we set 



0o = 0i = 0. (1.2) 
In the zeroth order we obtain a Vlasov-Poisson System 

S(fo) - Vz 02 • Vp/o = 0, /i = y fo dp, 
<j) 2 (t,x) = -f \z\~ 1 fJL (t,x + z)dz, f o (0,x,p) = f°(x,p) 

where 



(VP) 



S = d t + p ■ V x 

is the non-relativistic free-streaming operator. In the first order a linearized Vlasov-Poisson system 
appears 

S(fi) - V.03 • V p / - V,0 2 • V P h = 0, 

A*l = / /l rf P, A :r 3 = /i! + <9 t 2 0i. 



Hence, if we suppose that /i(0, x,p) = 0, we can set 

/i = and 03 = 0, 

which also yields 

/i! = 0. (1.3) 
In the second order we derive an inhomogeneous Vlasov equation coupled to a Poisson equation 

#(/ 2 ) - V x 02 ■ V p / 2 - V x 04 • Vp /o = 

2 2 

4/o5(0 2 ) + V -p ■ V,/ + <S(0 2 )p - V - V,0 2 ) - Vp/o I (LVP) 
/i 2 = y (/ 2 - l/2p 2 / ) dp, A x 4 = // 2 + «9 2 2 

for which we choose homogeneous initial data 

/ 2 (0,x,p)=0. (1.4) 

At this point we need to discuss the Poisson equation for 4 . Since A x df(f>2 = 47r<9 2 /io the Poisson 
equation in (|LVP|) can be rewritten as 

04 = A; 1 (47r/i 2 ) + (A x )- 2 (A7rd^ ). (1.5) 

Therefore, we define 04 by 

(f) 4 (t,x) = - J J \z\~ 1 (f 2 -l/2p 2 f )(t,x + z,p)dpdz-^ J \z\o$fo{t,x + z,p)dpdz. (1.6) 

Using (jVP|) and partial integration the second term can be rewritten as 

~ J \z\d 2 f (t,x + z,p)dpdz = ~ J J \z ■ p)d t fo(t,x + z,p) dpdz, z = z\z\~ 1 . 

Thus far this term is merely bounded instead of decaying at infinity and 04 seems to be deter- 
mined only modulo a function depending on t. But if the boundary condition is fixed in terms of 
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integrability in weighted Sobolev spaces (for a definition see section |2J below (|2.3|) ) instead of a 
pointwise estimates, the corresponding condition for vanishing at infinity is 

04 G L p s for every s < -2 + 3/p' (1.7) 

where 1 < p < oo, 1/p+l/p' = 1. We will show in section |2] that in fact 04 fulfills this integrability 
condition. 

The first aim of this paper is to show that 

f D := /o + c" 2 /2, (1.8) 
<P D = c- 2 (f) 2 + c- 4 (j) 4 

yields a higher order pointwise approximation of the Vlasov-Nordstrom than the Vlasov-Poisson 
system. We will call this system the Darwin system because in the case of individual charged 
particles the relevant approximation is usually called the Darwin approximation. It is clear that 
for achieving this improved approximation property also the initial data of the Vlasov-Nordstrom 
model has to be matched appropriately by the data for the Darwin system. For a prescribed initial 
density f° we are able to calculate (/o, 02) and (/2, 04) according to what has been outlined above. 
We then consider the Vlasov-Nordstrom system with the initial condition 

f(0,x,p) = f°(x,p) (1.9) 
(j>°(x) = (/>(0,x) = c~ 2 2 (O,a;) + c" 4 04(O,x) 
(j>\x) = d t <j>(0 1 x) = c- 2 ^0 2 (O,a;) + c- 4 ^04(O,a;). 

Before we formulate the main theorem of this paper let us recall that solutions of the Vlasov- 
Nordstrom system together with the initial conditions (jl.9|) exist at least on some time interval 
[0, T] which is independent of c > 1; see Thm. 3]. This time interval is fixed throughout this 
paper. 

Theorem 1.1 Assume that f° G C°°(1R 3 x R 3 ) is nonnegative and has compact support. From 
f° calculate (fo,4>2) an d (^2,04); and then define initial data for the Vlasov-Nordstrom system by 
M.y\) . Let (f, 0) denote the solution of the Vlasov-Nordstrom system with initial data M.ty) and 
let (f D , <p D ) be defined as in lil.ty) . Then there exists a constant M > 0, and also for every R > 
there is Mr > 0, such that 



\f{t,x,p) - 


-f D 


(t,x,p)\ 


< 


Mc~ A 


(x G M 3 ), 


\(j)(t,X 


)-< 


b D (t,x)\ 


< 


M R c~* 


(\x\ < R) 


\d t (f)(t,x) - 


-d t c 


b D {t,x)\ 


< 


Mc~ 4 


(x G M 3 ), 


\V x (f)(t,x) - 




b D (t,x)\ 


< 


M R c- 6 


(\x\ < R) 



for all p G R 3 , t G [0, T] and c > 1. 

The constants M and Mr are independent of c > 1 but do depend on the initial data. Note that 
if the Vlasov-Nordstrom system is compared to the Vlasov-Poisson system only, one obtains the 
estimates \f(t,x,p) - f (t,x,p)\ + c 2 \<p{t,x) - 2 (t,x)| + c 2 |V x 0(t,x) - V x 2 (t,x)| + \d t (j>(t,x)\ < 
Mc' 1 ; see Thm. 3]. 

While this approximation has the big advantage that, since by now the Vlasov-Poisson system 
is well understood, the existence of (/ , 02) and also of (/ 2 , 04) does not pose serious problems; 
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note that in (jLVPJ) the equation for / 2 is linear. Therefore one can hope to get more information 
on (jVNc|) by studying the approximate equations. 

As a drawback of the above hierarchy one has to deal with two densities fo, fi and two fields </>2 
and 04 to define f D . Thus it is natural to look for a model which can be written down using only 
one density and one field. The most natural candidate for such a system might be the following. 



*/+Pll-£)- 



^ = 



S(<P)p + c 2 (l 



2c 2 ) 



V p f = 4S(</>)f 



l -7^)f('r,P)dp 
An , Air . 2 _ 2 



Note, that (f D , 4> D ) solves this system up to an error of order c~ 4 . But since second derivatives 
occur it is not clear which initial conditions are to be posed. It turns out that it is more convenient 
to rewrite this system in terms of derivatives of <ft then in terms of 0. To this end we introduce the 
scalar force field ip* corresponding to dt<fi and the vector field E* corresponding to V0. This leads 
to the following system which we call Darwin Vlasov-Nordstrom system. 



An 



{^* + p-E*)p- 

7^)n,;P)dp, j* 



h) E \ 



c 2 AE* 



Pf*(-,-,p)dp 
47rV/i* + d t Vi/j* 



(DVNc) 



Again (f D , dt(f) D , V x (j> D ) solves ()DVNc|) up to an error of order c 4 

i 3 x 1R 3 

there exist c* > 1 and T* > such that the following holds for c> c 



Theorem 1.2 Assume that f° G C 00 ^ 3 x IR 3 ) is nonnegative and has compact support. Then 



(a) If there is a local solution of \D l/TVcj) . then the initial data ip° and E° of \D l/7Vcj) at t = 
are uniquely determined by the initial density f°. 

(b) The system \D VNc\) has a unique C 2 -solution (f*,i[)*,E*) on [0,T*] attaining that initial 
data (f°,ip°,E°) att = 0. 



(c) Define (ft* by 



^a-v-^a-(v-%*). 



;i.n) 



Then c 2 A0* = 4vr/i* + d t ip* and V x <f>* = E* as well as d t <p* = ip* + 0(c" 4 ). 

Let (f,<p) denote the solution of \ VNc\) with initial data (f°, 0*(O, •), d t <f)*(0, •)). Then there 
exists a constant M > 0, and also for every R > there is Mr > 0, such that 

\f(t,x,p)-f*(t,x,p)\ < Mc~ 4 (x G R 3 ), 

\V x (f>{t,x)-E*(t,x)\ < M R c~ G (\x\<R), 

\d t (j)(t,x) -i/;*(x,t)\ < Mc~ 4 (xGM 3 ), 

\<j>(t,x)-<p*(t,x)\ < Mc~ A (xeR 3 ), 



for allpe R 3 , t G [0, min{T, T*}], and c>c*. 
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We shall continue with a few remarks about the Darwin approximation of the relativistic Vlasov- 
Maxwell system derived in [2|. In that paper only the 1 PN approximation is treated but we want 
to emphasize that in the case of only one species it is possible to establish the 1.5 PN approximation 
for Vlasov-Maxwell as well in the same way as it is done for (jVNc)) in this paper. On the other 
hand, if several species with different charge to mass ratio, say / + with charge 1 and /~ with 
charge -1, are under consideration, mass is set to unity for both species, an effect due to radiation 
damping takes place, which is proportional to 



-3i 




dt VU Q (f+ + fo){t,x,p)dpdx (1.12) 



(VPplasma) 



where (f^ 1 , Uq) is a solution of the Vlasov-Poisson system in the plasma case with two different 
species 

± V X U ■ V p / ± = 
AU = 4n [(ti-fo)dp. 

(In the case of only one species and in the gravitational case the corresponding term vanishes, see 
lETBI) below.) Note that d t J J Vf/ (/ + + fo)(t,x,p) dpdx = ^D{t) where 



D(t)-= x(f+-fu)(t,x,p)dpdx (1.13) 




is the dipole moment. Thus it seems reasonable to replace ( VPplasma ) by 

d 3 



± (^D(t) + v x u ) ■ v p /± = 

AU = An J(ti-fo)dp. 
D(t)= / x(/ + - fo)(t,x,p)dpdx 



(VPrad) 



where a ^ is a constant which has to be calculated from the full system. However, in this system 
phase space has increased by 3 dimensions and it is not so clear how to single out the initial data 
of -^D, note that D(0) and J|-D(0) are already determined by the initial data of f ± . On a formal 
level (jVPrad|) can be reduced in several ways, which are discussed in and [TT] . 

In the gravitational case effects due to damping are expected to take place in the 2.5 PN 
approximation and to be connected with the fifth time derivative of the quadrupole moment, see 
|llj . Thus, it seems reasonable that approximations up to the order of 2PN conserve an energy 
which can be calculated by an expansion of the energy 



£ = c 2 J J Vl + c- 2 p*f(t,x,p)dpdx + c 2 J (\d t <f){t,x)\ 2 + c 2 \V x <j>{t,x)\ 2 } dx, (1.14) 



which is conserved by solutions of the full Nordstrom system; note that in the Vlasov-Maxwell case 
the corresponding statement is true for approximations up to the 1PN level. However, here this 
fails already in the Newtonian approximation in a rather obscure manner; the conserved energy 
of (jVPj) is J J y/ dpdx — J \ V x (p2\ 2 dx whereas from the expansion of ()1.14j) one would expect 
J J y/ dpdx + J \V x 4>2\ 2 dx to be conserved. But this is the energy of Vlasov-Poisson system in 
the plasma case. 
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The paper is organized as follows. Some facts concerning (|VPJ) . (jLVPJ) . and (jVNc|) are collected 
in Section |21 as well as the proof of (jl.7j) . The proof of Theorem 11.11 is elaborated in Section El 
Section 0] contains the proof of Theorem 11.21 For these proofs we will mostly rely on suitable 
representation formulas for the fields (refined versions of those used in jHO]), which are derived in 
the appendix, Sectional 

Notation: B(0,R) denotes the closed ball in R 3 with center at x = or p = and radius R > 0. 
The usual L°°-norm of a function = 4>{x) over x G R 3 is written as \\4>\\ x , and if = <f>(x,p), we 
modify this to H^H™. For m G N the W m '°° -norms are denoted by ||0|| ma ., etc. If T > is fixed, 
then we write 

g{t,x,p,c) = O cpt (c- m ), 
if for all R > there is a constant M = Mr > such that 

\g{t,x,p,c)\ < Mc' m (1.15) 

for |x| < R, p el 3 , t G [0,T], and c > 1. Similarly, we write 

g(t } x,p,c) = 0(c- m ) 

if there is a constant M > such that (|1.15|) holds for all x,p G M 3 , t G [0,T] and c > 1. In 
general, generic constants are denoted by M. 

2 Some properties of (EE), JEEEDj and dYHH) 

There is a vast literature on (|VP|) . see e.g. [7| Sect. 4] or [21] and the references therein. For our 
purposes we collect a few well known facts about classical solutions of (|VP|) . 

Proposition 2.1 Assume that f° G C 00 ^ 3 x R 3 ) is nonnegative and has compact support. Then 
there exists a unique global C 1 -solution fo of I VP}) , and there are nondecreasing continuous func- 
tions Pvpi Kyp : [0, oof— ► R such that 

\\fo(t)\\ x , P < \\f%, P , 

supp/ (t, v ) C 5(0,Py P (t)) X S(0,Pyp(t)), (2.1) 

ll/o(*)lli^p+ ll^(*)Ila^ < ^vp(*), 

for t G [0, oof. 

This result was first established by Pfaffelmoser [20 J, and simplified versions of the proof were 
obtained by Schaeffer [21] and Horst |H] ; a proof along different lines is due to Lions and Perthame 

For our approximation scheme we also need bounds on higher derivatives of the solution. This 
point was elaborated in JHI where it was shown that if f° G C fc (R 3 x R 3 ), then / possesses 
continuous partial derivatives w.r.t. x and p up to order k. The existence of continuous time- 
derivatives then follows from the Vlasov equation. Thus, / and 2 are C°° if f° is C°°, and by a 
redefinition of Kyp we can assume that 

II/o(*)|| w <^vp(*), tG[0,oo[. (2.2) 
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The existence of a unique C 1 -solution f 2 of (|LVP|) follows by a contraction argument, but we 
omit the details. Furthermore it can be shown that there are nondecreasing continuous functions 
Plvp-i Klvp '■ [0, oof— > R such that 

Bupp/ 2 ( v ,t) C B(0,P LVP (t)) x B(0,P LVP (t)), (2.3) 
ll/ a (*)ll af . lP + W^^^vp(*). ( 2 - 4 ) 

for t G [0, oo[. 

Next we shall investigate the boundary behavior of 04 in terms of integrability in weighted 
Sobolev spaces where the weight function is defined by p(x) = (l + lx] 2 ) 1 ^ 2 . For G No, 1 < p < oo 
and s G R we say u G W*> p if and only if u G W£f(R 3 ) and p s+ l Q l9 a M G L P (R 3 ) for < |a| < fc; 
a G Nq a usual multi-index. In order to establish (J1.7)) we use some mapping properties of the 
Laplacian in weighted Sobolev spaces due to McOwen. Citing [TH], for 1 < p < oo, l/p+ 1/p' — 1, 
we have that 04 G for s < — 2 + 3/p' if <9 t 2 2 G £ 2+s . Furthermore, <9 2 2 G L 2+s holds true if 

A<9 2 2 = d 2 fx G L^ +s and ^ h(x)d 2 fi (x) dx = (2.5) 

for every polynomial ft of degree less than 1. Thus, to prove (jl.7|) it is sufficient to show (|2.5|) . 
Since <9 t 2 /io(£, ■) has compact support for every £ > 0, compare (j2.1j) . the first condition is clear. 
Now let ft(a;) = a + aiXi + a 2 x 2 + a 3 x 3 , a.j G R. Employing (jVP|) . some partial integrations yield 

h(x)d 2 fi (t,x)dx = —J J h(x)p ■ V x d t fo(t,x,p) dpdx 

= Vh- J J pd t f (t,x,p) dpdx = Vh ■ J J p(-p ■ V x f + V x fa ■ V p / )(£, x,p) dpdx 

= — V/i ■ J V(f)2^o{t,x) dx = — — a,i J difodjdjfoit, x) dx 

= — V/i • [ V\V(p2(x)\ 2 dx = 0. 
8tt y 

Concerning solutions of (jVNc|) we have the following from j3J Thm. 3, Proposition 2]. 

Proposition 2.2 Assume that f° G C°°(R 3 x R 3 ) is nonnegative and has compact support. 7/0° 
and 1 are defined by il.ty) . then there exits T > (independent of c) such that for all c > 1 the 
system h VNc}) with initial data hi. 9^1 has a unique C 1 -solution (/, 0) on £/ie time interval [0, T]. In 
addition, there are nondecreasing continuous functions (independent of c) Pvn,Kvn '■ [0,T] — > R 

/(z,p,t) = if \p\>P V N(t), (2.6) 
|/(t)U P +|0(t)| M <tf w (*), (2.7) 



for all x G R 3 , t G [0,T], and c > 1. 



3 Proof of Theorem 11.11 

The proof follows the usual lines developed in (22] and in [2| for higher order approximations. In 
a first step the difference of the fields, — <fP and the difference of its derivatives is estimated in 



S 



terms of the difference of the densities h = f — f D . For that purpose we shall use quite elaborated 
representations of the fields, which are derived in the appendix. Next we will calculate which Vlasov 
equation it is that h fulfills. Using the matching of the initial data and a Gronwall argument, it 
follows that \h\ = 0(c~ 4 ) which in turn yields the announced estimates of the errors in the fields. 

As V0 enters the Vlasov equation with the factor c 2 we have to be most precise in comparing 
V0 with V0 D . Therefore, we give the reasoning for that term in some detail. In section 15.1.11 
below we will show that the gradient of the approximate field (f> D from (|1.8|) admits the following 
representation. 

V,0 D = < ext + < int + < bd + O cpt (c-') (3.1) 

with 



+ 2? 

1 
1 

+ - 

c 3 

1 



0*,bd(^) = -77 



1 

+ ? 

1 

1 

cH 

t 
3c* 
1 
3c^ 



\z\>ct 

r 

\z\>ct 
z\<ct 
z\<ct 
z\<ct 
z\<ct 
z\<ct 
z\<ct 



\A 2 z(j D - ^f^{t,x + z,p)dpdz 
zd 2 fo(t, x + z,p) dp dz, 



(3.2) 



- 2 zf D (t(z),x + z) dz 



-2 



(2z(z ■ p) — p) f D (i(z), x + z,p) dpdz 



~ 2 ( - 3z(z ■ p) 2 + (z ■ p)p+ -zp 2 ) f D (i(z),x + z,p) dpdz 
~ l zz ■ V x 4>2f D (t( z )i x + z iP) dpdz 



A(z ■ p) 3 z — (z ■ p) 2 p — 4(z • p)p 2 z + p 2 p)f v (t(z),x + z,p) dpdz 



-1; 



z[2(z- p)z ■ V x fa - p ■ V x <j> 2 - S(<p 2 ))f D (t(z),x + z,p) dpdz 



z <ct 



dt[ V x <p 2 ^o)(i(z),x + z) dz, 



(3.3) 



z\=ct 



(z ■ p) (z • p) 2 — p* 



z \=ct 



z\=ct 



z(z-p)[l 
z(z ■ p)d 2 fo(0, x + z,p) dpds(z) 
(z ■ p)pd t f (0, x + z,p) dpds(z) 



f°(x + z,p) dpds(z) 



where the subscripts 'ext', 'int' and 'bd' refer to the exterior, interior and boundary integration 
in z. We also recall that z = \z\~ x z and t(z) = t — c~ 1 \z\ is the retarded time. On the other hand, 
according to Section 15.1.21 below, we have 



x,ext 



+ 



+ O cpt (c 



(3.4) 



with 



Iff/' t 2 t 3 \ 

4x,*tt(t,x) = — 2/ / \z\~ 2 z [fo + td t f + -d 2 f + -dff )(0,x + z,p) dpdz 
c J\z\>ctJ v 2 6/ 
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+ 2^i / / \z\' 2 zp 2 (f + td t f )(0,x + z,p)dpdz 
+^J n>ct J z($fo+t$f o yo,x + z,p)d P dz 
~ 2 zd t j 2(0, x + dpdz, 



+ - 



+ - 



1 
1 



>rf 



<rf 



<ct 



<ct 



(3.5) 



-2- 



zf(t(z),x + 

~ 2 (2z(z • p) — p)f(i(z),x + z,p) dpdz 

3 

~ 2 ( - 3z(z • p) 2 + (z ■ p)p + - zp 2 ) f(t(z),x + z, p) dp dz 
~ x zz ■ V x (j)f(t(z),x + z,p) dpdz 

(a(z ■ p) 3 z — (z ■ p) 2 p — 4(z ■ p)p 2 z + p 2 p^ f (i(z),x + z,p) dpdz 



-2 



1 z\2(z -p)z ■ V x (j)-p - V x (j) - S{4>)jf{i(z),x + z,p) dpdz 
-3^5 J < d t (v x (P 2 ^)(0,x + z)dz, (3.6) 



(j>x,m(t,x) = (j) x>bd (t,x). 

In order to verify (jl.lO|) we start by comparing the exterior fields. Let x £ B(0, R) with R > be 
fixed. Then we obtain from ()3.2|) and ()3.5j) . taking into account that 



f 2 (0,x,p) dp = 



by (frVPf and dUJ) as well as fl2~Tj) . (Q , (Q and (JH, 



|V x ai t(t,a;) - V a .0 6xt (t,x 

-2 



< 4 / l~l 

C -'|z|>Ct 



M 
< 



t 2 t 3 

fio(t, x + z) - (iiq + £<9 t // + — d 2 fi + —dffio) (0, a; + 



d 2 /j>o(t, x + z) — (d 2 fXo + tdfuo) (0, x + z) dz 
+ \ I I \z\~ 2 \f 2 (t,x + z,p) - (f 2 + td t f 2 )(0,x + z,p)\dpdz 

c J\z\>ctJ 

I I \ z \~ 2 P*\fo(t> x + z >P) ~ (/o + td t f o )(0,x + z,p)\dpdz 

ZC J\z\>ctJ 

(J (t - s) 3 P V p(s) 3 K V p(s)l Bi0 ,p V p(s))(x + z) ds^j dz 
+ ~T / ( / (* _ s)P V p(s) 3 K V p(s)l B(0t p vp{s)) (x + z) ds ] dz 

c J\z\>ct \J0 J 



dz 
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M 



M 



/ 


z\ 


l\z\>ct 




f 


z\ 


l\z\>ct 




1 


z\ 


>\z\>ct 





J (t - s)P V p(s) 5 K V p(s)1b(o,p vp ( s ))(x + z) ds*j dz 

b(o,p lvp (s))(x + z) ds^j dz 

Mt 4 f Mt 2 f 

< — 2~ / \z\~ 2 1b(o,r+m q )(z) dz ^ — (\z\~ 2 + l)l B(0}R+Mo )(z)dz 

c J\z\>ct c J\z\>ct 

< M R c- 6 ; (3.7) 

note that here we have used 

M = max (p V p(s) + K VP (s) + P L vp(s) + K LVP (s)) < oo 
se[o,T] \ J 

and for instance 



„2 



/ \z\ 2 1b(o,r+m ){x + z)dz < ! \z\ 2 dz < M R c 

J\z\>ct c t J\z\<R+M 



To bound \<f> Xj i n t{t, x ) ~ 0fint(^ a; )l we shall first treat the last term in (|3.3|) and (|3.6|) respec- 
tively, which is special in some way as it does not arise from the expansion of the integral kernels 
of V x <p. Regarding this term, a crucial observation is that by (|VP|) and partial integration 



V x 2 M^ x ) dx = j~ ( S / 9 ^d Xj d x .(f)2{t, x) dxj 

V - =1 J J i=l,2,3 

= d Xi d x .(j) 2 d x . ( p 2 (t,x)dx^ = ~J V x \V x M^ x )\ 2 dx = (3.8) 

for all t G [0, oof. This yields 

_ 3~5 / d t ^ V x (f)2^(t(z),x + z)dz = -— -g. I d t ^fi V x (p2^(t,x + z)dz + O cp t(c' 6 ) 



" <TF (/ {^(Nx4>2j (t, x + z) dzj + 3-4 / (aioV x 2 j ( 0. : ) dM : ) £>,.,,/ ( r" 



6> 



= A( / {j^x<h){t,x + z)dz\+^- i f (jL V x <l)^){Q,x + z)ds{z) + O cpt (c- 6 ) 
= ^s/ > dt(lM)Vx<h)(t,x + z) dz + O cpt {c- & ) (3.9) 
where we used that for x G 5(0, i?) 

/ dJ/ioVzcM^^x + z) -<9 4 (/i O V a; 02)(^a; + ^)^ 

J|z|<rf v 7 v 7 

< — I \z\K^ P (t)P^ P (t)l B ( ,p VP )(x + z) dz 

C J\z\<ct 



M 

< — 

c 



J \z\l R+Mo {z) dz < M R c 1 
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(note that \t(z)\ < t for \z\ < ct and K V p and Pyp are nondecreasing) . Again by (J3.8|) and same 
calculations we also have 

~3c5 / fy(j^V x fa)(Q,x + z)dz = -^ J d t (jiq V x 2 ) (0, x + z) dz. (3.10) 
Employing ()3.9|) and ()3.1U|) we can estimate 

d t (j"oV^ 2 ) {i{z),x + z) - d t \^oV x (j)2j (0, x + z) dz 

J | d t (jM)V x (h) (t, x + z)-d t (jioVxfa) (0, x + z) dz + M R c- 6 

[ \z\K 2 VP (t)P VP (tfl B{0)Pvp{t) (x + z)dz + M R c~ 6 < M R c~ 6 (3.11) 

J\z\>ct 



1 

3c5 



z\<ct 





M 


< 






IF 




M 


< 









for all t E [0, T] and x E 5(0, 

Next we recall from [3J Thm.3] that 



|V^(t,z) - c~ 2 V x (j>2(t,x)\ = 0(c 



(3.12) 



actually in [3] the initial conditions are different, but we only added terms of order c -4 so that an 
inspection of the proof in [3] leads to (|3.12|) . Furthermore, in section 15. 1.11 ()5.5|) and section 15. 1.21 
(|5.31|) it is shown that 

\d t( j> D (t,x)\ + \d t ^t,x)\=0(c- 2 ). 



Also, we define 

H(t) = sup{\f{s,x,p) - f (s,x,p) - c- 2 f 2 (s,x,p)\ : x E IR 3 , p E IR 3 , s E [0,t]}. 
Now we have to proceed in two steps: Using (|3.12j) and (|3.13|) prove 

|V x D (t,x)- V,0(t,x)| < \M R (c- 3 + H(t)) 

c z 

\dt<p D (t,x)-dt<j)(t,x)\ < M(c~ 3 + H(t)), 
H(t) < Mc~ 3 



(3.13) 
(3.14) 



employ this to derive 
for all t E [0, T], which in turn gives 



\V x <t> D {t,x)- V x (f>{t,x)\ < M R c~ 5 
\d t (f) D (t,x) - d t (f)(t,x)\ < Mc~ 3 



and especially as c 4 V X 02 = 0(c 



-4\ 



\c~ 2 V x (j) Q {t,x) - V x <f>{t,x)\ < M R c~\ 



(3.15) 



In the second step this procedure will be repeated where estimate ()3.12j) is replaced by 1)3. 15J) . 
which will give the announced estimates. To avoid redundancies we shall only carry out the second 
step, and therefore assume that ()3.15|) is already proved. 
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Introducing the constant 

Mi = max ( P VN (s) + P V p(s) + Plvp(s)) < oo 

se[o,T] V ' 

it follows that f(s,x,p) = f (s,x,p) = f 2 (s,x,p) = for x G M 3 , \p\ > Mi and s G [0, T]. Also if 
i?o > is chosen such that f°(x,p) = for |x| > Rq, defining the constant 



M 2 = Rq + TM 1 + max (p vp {s) + P LV p(s)) < oo, 
se[o,T] V / 



(3.16) 



we have that f(s,x,p) = f (s,x,p) = f 2 (s,x,p) = for |x| > M 2 , p G M 3 and s G [0,T]. Let 
again x G B(0,fl) with i2 > be fixed. From (f3~HJl . (f3~TTJl . (pTSjl . (j3~T2j) . (jHHSj) (fOjl and 

< t(-z) < t for |z| < ct we obtain 



: (t,x)-0^ nt (t,x)| 



1 



1 

+ - 

& 

1 

+? 

1 

+ ^7 



1 



|z|<ct 

r 

\z\<ct 

f 

\z\<ct 

f 

\z\<ct 

f 

\z\<ct 

f 

\z\<ct 



-2 



(f-f°)(t(z),x + z,p) 



dpdz 



\2z(z-p)-p\Uf-f D )(t(z),x + z,p) 



dpdz 



3z(z-p) +(z-p)p+ -zp | (/ - / )(t(z),x + z,p) 



dpdz 



(V a 



f(t(z),x + z,p) dpdz 



v x Mf-r)(t(z),x + z, P )) 



dpdz 



M 

+— 

c d 

M 

+ 7 



\z\<ct 

f 

\z\<ct 



-2 



i-l 



3 p 2 

A{z -pfz + (z-p) 2 p - A(z-p)p 2 z + — p (f - f D )(i(z),x + z,p) 



dpdz 



\P\ 



V x 5-V x 02 



+ 



S(<p) - \s(fc) )f{t{z),x + z,p)dpdz 



D 



(t(z), x + z,p) dpdz + M R c 



< ^(M? + M?)H(t) [ \z\- 2 l B{0M2) (x + z)dz 

C J\z\<ct 

+^(M 3 + M$)K VP {T)H{t) I W^b^m^x + z) dz 

C J\z\<ct 



M, 



^{Mt + Mt)\\f°\\ x , p I l B{0>M2) (x+z)dz + 



M, 



R 



\z\<ct 



< c- 2 M R (H(t) + c- 4 ) 
since for instance 



(3.17) 



/ \A 1 l B (fi,M2)( x + z ) dz < / \A 1 dz<M R 

J\z\<ct J\z\<R+M 2 



\z\<R+M 2 

Recalling that Xjb d = 0^bd> we can summarize (|3.1|) . (|3.4|) . (|3.7|) and (|3.17|) as 

c 2 \V x <P(t,x) - V x <p D (t,x)\ < M R (H(t) + c" 4 ) 



(3.18) 
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for x G 5(0, R) and t G [0, T}. Formulas (I5U). (ET2%J) . (ETHJ) . (ICTIj) and an analog calculation also 
lead to 



\(f){t,x) -<p D (t,x)\ < M(H(t) + c- 4 ) 
\dt<f>(t,x) -d t (f) D (t,x)\ < M(H(t) + c- 4 ) 



(3.19) 
(3.20) 



for t G [0, T]. It remains to estimate h = f- f D . Using (jVNcjl . Q|) . dVP) and (TCVPl . it is found 
that 



5(0)p + 7 c 2 V^ 



-5(/ D ) + 5(0)p + 7 c 2 V x -V p r + 45(0)/ 



1 2 ' 

P - P + -j^V V 



+ s(<f>)p-V p f 2 + 



V <B /o-^(p-p)-V !B / 2 + 5(0) 



1 5(02)1 p ■ Vp/o 



7C 2 V X (; 



2c 2 / 



+ 



7Vx0 - ^V x 2 



V x .0 2 - -V 

1 X 



V p / 



V p / 2 + 45(0)/i - 4 (-5(0 2 ) - 5(0)) /o + ^5(0) / 2 . 



If |p| < Mi, then also |p| = (1 + c 2 p 2 ) 1//2 |p| < |p| < M\ uniformly in c, and hence 



p — p 



p-p+ ^P 2 P 



< Ma- 



in view of the bounds fl2UJ, and fl27ZJ), 5(0) = C(c~ 2 ) as shown in section EH (1q~23J) and 
(JOIl), thus by (|3~T%1) and (l3~27H) 



S(h)(t,x,p)- S((f))(t,x,p) +7c 2 V x 0(t, x) -V p f(t,x,p) 
< M( C - A + H{t)) 



(3.21) 



for all |x| < M 2 , |p| < Mi and t G [0,T]. But in {(t,z,p) : \x\ > M 2 } U {(t,a;,p) : \p\ > Mi} we 
have h = f — f D = by the above definition of Mi and M2. Accordingly, ()3.21|) is satisfied for 
all x G M 3 , p G IR 3 and t G [0,T]. Since fi,(0,x,p) = the argument from [3| (see also [231 p.416]) 
yields 



H(t) < [ M{c~ 4 + H{s))ds, 
Jo 



and therefore H(t) < Mc~ A for t G [0,7]. Then, due to ff3~T3j) - ([3~2Uj ) c 2 |V x 0(t,x) - V z D (t,x)| < 
M^c" 4 for |ac| < -R and t G [0, T] as well as |0(t, ac) - D (£, ar)| + \d t (j)(t, x) - d t <p D {t, x)\ < Mc~ A for 
x G IR 3 and t G [0, T\. This completes the proof of Theorem ll.il □ 



4 Proof of Theorem 11.21 

In this section we will be sketchy and omit many details, since the proof is more or less a repetition 
of what has been said before. First let us assume that there is a C 2 -solution (/*, ip*, E*) of (|DVNc|) 
existing on a time interval [0, T*} for some T* > such that supp f*(t, •, •) C IR 3 x IR 3 is compact 
for all t G [0,T*]. Then, after partial integration 

tp*(t,x) = — I I \z\~ 2 z ■ pf*(t,x + p,p) dpdz (4.1) 
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Air 

AE*(t,x) = — (V/i* + $W). 



(4.2) 



C 2 47T 
1 

1 

1 

2^ 



Since f*(0,x,p) = f°(x,p), it follows that ip*(0, •) is determined by f°. 

In order to compute the Poisson integral for E*, we calculate by means of the transformation 
y = w — z, dy = dw, and using ()5.35|) below, 

\x-y\- l Vd t r{t,y) dy 

\ x -y\~ l ^y(^J J \ z \' 2 ( z -P)9tf*(t,y + z,p)dpdz^j dy 
\z\~ 2 z\x — w — z^ 1 dz^j J pd m d t f*(t,w,p) dpdw 
\w — x\^ 1 (w — x) ■ pV w d t f*(t, w,p) dpdw 

2c 1 / J( z 'P) V ^f*( t ^ x + z ^P) d P dz 

1 ( — z(z ■ p) + p)d t f*(t, x + z,p) dpdz. 

If we invoke the Vlasov equation for /* and integrate by parts, this can be rewritten as 



i=l,2,3 




1 

2? 




1 

2c 1 

+ 




z(z-p)+p)[ -(l-^)p.VJ 



(r+P- E*)p + c 2 (l - 7^)E*] ■ V P f* + 4(^* + p ■ E*)f*j (t, x + z,p) dpdz 



2c 1 / /'^' 2 (~ 2 ( z 'P">P + 3 ( z 'P") 2z ~P 2z )( 1 ~ 7^)f*( t ' x + z 'P) d P dz 
V-z®z) [O* + p ■ E*)p + c 2 (l- |^) £*] f* 

+ [z(z ■ p) — p) (ip* + p ■ E*) f* | (t, x + z, p) dp dz. 



2c 4 




Therefore the solution E* of ()4.2|) has the representation 

E*(t,x) = A-^^V/ij+A- 1 ^^) 
= — - J \z\~ 2 z/j*(t, x + z) dz 

+ h S I \ z \~ 2 {- 2 ^-p)p +? >(z-p) 2 z-p 2 z) (i-^jr{t,x+z,p)d P dz 



1 

'2c 1 




z\- L Ul-z®z) (^* + p ■ E*)p + c 2 1 



PL\ 
2c 2 ) 



E* 



r 
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+ [z{z ■ p) — p) (ip* + p ■ E*) f* | (t, x + z, p) dp dz. 



(4.3) 



In particular, if we evaluate this relation at t = the Banach fix point theorem applied in C&(R 3 ) 
shows that for c > c* sufficiently large the function E*(0, •) is uniquely determined by f° = 
/*(0, •)■ Thus f° alone fixes ^° and E° . 

Concerning the local and uniform (in c) existence of a solution to (jDVNc|) one can use (J4.1)) 
and (|4.3|) to follow the usual method by setting up an iteration scheme for which convergence can 
be verified on a small time interval: cf. [7J Sect. 5.8] and also 4J. For comparison of E* with V x (j) D 
we give an alternative expression of V x <f) D . Using (jVP|) and partial integration twice we have 




-V x ^j J \z\dtfi (t,x + z) dz = ^J zd 2 /j, (t,x + z) dz 

zd t ^-p- VJo + V x .0 2 ■ Vp/o) {t, x + z,p) dp dz 
- (z ■ p)z + pjd t fo(t,x + z,p) dpdz 
l (-{z-p)z+p\(-p- V x f + V x 4> 2 ■ Vp/oJ (*, x + z, p) dp dz 
2 (y — 2(z • p)p + 3(z • p) 2 z — p 2 z^j f (t, x + z,p) dpdz 
| ^ | ~ x ( 1 - z ® z)V 'vfafiofax + z) dz. (4.4) 





1 

2^ 
1 

2^ 
1 

2^ 
1 

2c* 

1 

Together (JBTTJl . (IQ1) (jOJ) and (|5^jl . (jUJ) reveals the analogy of to V X D and ^* to <9 t D 
respectively at the relevant orders of c _1 . Comparison of with ()VP|) . (jLVP|) and ()1.8|) gives 

the analogy of the initial values of (jVNcjl . Finally, by similar arguments as used in the proof of 
Theorem II. II it can be shown that solutions of (jDVNcJ) approximate solutions of (|VNcJ) up to an 
error of order c~ 4 . 




5 Appendix 

5.1 Representation Formulas 

5.1.1 Representation of the approximation force field 

Here we will present representation formulas for the approximate field and its derivatives. As the 
calculations are lengthy and arguments are similar we shall only sketch the computations leading 
to the formula for V X 4> D and merely give the formulas for (fi D and dt<j) D . 
We recall (fi D = c~ 2 (p2 + c _4 04 where 

<j) 2 {t,x) = — I \z\~ l ji Q {t,x -\- z) dz 



4>i(t,x) = — — J \z\d 2 no(t, x + z) dz — J \z\ 1 n 2 (t, x + z) dz, 
thus, using //„ = / /„(-, -,p) dp and // 2 = / (/ 2 - y/oj (•, -,p) dp 

V^ D (t,x) = J\ z \-*z(f D -^f o yt,x + z,p)dpdz 
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J zd%/j (t,x + z)dz. (5.1) 

We split the domain of integration in {\z\ > ct\} and {\z\ < ct}; note that the exterior part gives 
<t>x,ext- To handle the interior part {\z\ < ct} we expand the densities w.r.t. t about the retarded 
time 

t(z) := t - c~ l \z\. 

To begin with we have 



1 



\z\<ct 



z\ 2 zfi (t,x + z) dz 

~2 S 1^(1 + ^ + ^ + ^)f, (i(z),X + Z)dz 



\z\<ct 



— 7-0 / \z\ 2 z / (t — s) 3 dffi (s, x + z) ds dz (5.2) 

6c J\ z \<ct Ji(z) 

Since c\/z + / p • V x fo dp = by ([VPjh we also find 

/ \z\~ x zd t yLo{i{z),x + z) dz = / kT^f / P ■V x fo(i( y z),x + z,p) dpdz 

C J\z\<ct C J\z\<ct 

k I /M-V (v,W ( », I+ ,ri] + c-^/ D («( z) , I+ , P ))*^ 

- J\z\<ctJ V 7 

I (z • p)f°(x + z,p) dpds(z) 

- 2z(z ■ p) +p)f (i(z),x + z,p) dpdz 
z\~ l z(z ■ p) dtfo{t(z), x + z,p) dpdz, (5.3) 



c 
I 

~cH 

1 



|z|=et 



2 <ct 



|*|<cf 



observe that = for |z| = ct was used for the boundary term. Using (|2.ip . ()2.2|) and ()2.4|) the 
last term in ()5.3|) as well as the remaining terms coming from the interior part of (|5.1|) and from 
(|5.2|) are in O cpt (c~ 4 ); note that |x| < R for some R > together with the support properties of 
/o imply that we only have to integrate in z over a set which is uniformly bounded in c > 1. Thus 
the leading terms in the orders c~ 2 and c~ 3 are specified. The following scheme emerges, 

• expand the several terms in the interior part of (J5.1|) up to order c~ 5 , the error term is in 
O cp t{c 6 ); 

• replace time derivatives on the densities f and / 2 via Vlasov's equations by space derivatives 
according to ijVFj) and (jEVFJ); 

• carry out partial integrations, using 

(V x g)(i(z),x + z) = V z [g{Kz),x + z)}+c- l zd t g{i{z),x + z) 
in the case of derivatives with respect to x. 
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Following this scheme one can calculate one after another the leading terms up to order c 5 . 

Addition of some additional terms in C(c~ 6 ), as e.g. J\ z \< ct I \ Z \~ 2 ^P 2 h{^i z )i x + z,p) dpdz, 
then gives flUS)) and (JH3J). 



One can also calculate 
with 



+ & + <& + O(o' 4 ) (5.4) 



0?xt(^) = — 2 / N H^o + c 2 fi 2 )(x + z,t)dz 

C J\z\>ct 



1 

'2^ 



z[>ct 

|z|<9 t 2 /zo(i, a; + z) dz, 

|z|>ct 



0f nt (t,x) = -i / / |z|-7"(t>),x + z,p)dpdz, 



z|<et 
|z|=ct 



= +4 / I (z -p) f°(x + z,p)dpds{z) 



and 



^ = Cxt + Cnt + + 0(O (5.5) 
is derived in the same manner with 



<t>t,ext{t, X ) = ~2 / I \z\ 2 {z-p)f D {t, C+ Z.j>)(lp<lz. 



0fint(M) = 4 / / M 2 (z ■ p)f v (t( y z),x + z,p) dpdz 



c 

1 



c 3 



|z|>et 

r 

|z]<ct 

z|<Ct 



z| 1 z -V x 4)2f D {i{z),x + z,p) dpdz 
\ J f \z\~ 2 (p 2 -2(z.p) 2 )f D (i(z),x + z,p), 



'\z\<ct ■ 

= — it/ I (z -p) 2 f (0,x + z,p)dpds(z). 



cH 



\z\=ct 



5.1.2 Representation of the Nordstrom force-field 

We recall the following representation from |3J Proposition 3; Proposition 4] 

4> = 4>d + 4>s (5.6) 

9 t (j) = 4>t,D + 4>t,BD + 4>t,a + 4>t,b + 4>t,c (5.7) 

V x = <f) XtD + (j) Xt BD + 4>x,a + 4>x,b + 4>x,c (5.8) 



where 

(f) D {t,x) 



9t( — 4> Q {x + ctuj) dw ) + — / 4> 1 (x + ctuj) duj 
V 47r J\u,\=\ / 4vr J m=1 



z|<ct 

<Pt,D(t,x) = d t (p D (t,x) 



4>s{t,x) = — n / / M l{p)f{t(z),x + z,p)dpdz 
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and 



h,BD\ 



(t,x) = -c- l (ct)- 1 I bd <t,t f°( y x + z,p)dpds(z 



z\=ct 



<PtJt,x) = -c 2 / / a*f(i(z),x + z,p)dp °- 



2 

dz 
\z\ 

<>,Jt.-r) = r 1 / / C**(V x (j>)f(i(z),X + Z,p)dp^- 



'\z\<ct 

<P tfi (t,x) = -c- 2 [ I b^S(<f>)f(t{z),x + z,p)dp 

J\z\<ct 

f 

\z\<ct 



<t>x,D{t,x) 


= V,( 


i>D{t,x) 


>x,BD(t,x) 


= — c 


-Hay 


<j) x ,a(t,x) 


= —c 


in 

J\z\<ct 


4>x,b(t,x) 


= —c 


J\z\<ct 


4) x ,c{t,x) 


= —c 


i 

J\z\<ct 



where the kernels are 



bd^f°{x + z 1 p)dpds{z) 

\z\=ct J 

J a^ x f(i(z),x + z J p)dp-^ I 

dz 

b 4> *S(<i))f(i(z),x + z,p)dp — 

dz 

c^(V x (f))f(i(z),x + z,p)dp — 

z 



bd^ = 7 (p)(l + c~ 1 z-p)- 1 

= -^{l + c~ l z-p)~ 2 p-{z + c~ l p) 

b^ = -fil + c^z-p^iz + c^p) 2 

c 4>t = 7 3^ + c ~ l z-p)- 2 (z + c- 1 p) 



and 

M** = 'y(p)(l + c~ 1 z-py 1 z 

a <t>x _ ry(i _|_ c -ig . p) _2 |^ -\- c ~ l p — c~ 2 p A (zAp)} 
c ^ = z©c^GM 3x3 . 

Next we expand these kernels in powers of c™ 1 (in this subsection we will explain the argument 
leading to a formula for V x (j) in detail and only give the formulas for the other fields). According to 
we can assume that the p-support of f(t, x, ■) is uniformly bounded in i 6 R 3 and t G [0, T], 
say f(t,x,p) = for \p\ > P. Thus we may suppose that \p\ < P in each of the p-integrals, and 
hence also \p\ = 7(p)|p| < \p\ < P uniformly in c. It follows that 

p=(l-l/2c- 2 p 2 )p + 0(c' 4 ). 

For instance, for the kernel bd^ x of (f) Xi BD this yields 

bd** = ^7(p)(i + ^r 1 
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1 - l/2c~V + C(c~ 4 ) 1 - c~\z ■ p) + c~ 2 (z ■ p) 2 + c" 3 (l/2(z ■ p)p 2 - (z ■ pf^j + 0{c 
= z(l-c-\z-p) + c~ 2 ((z-p) 2 -l/2p 2 )+c~ 3 ((z-p)p 2 -(z-pf)^+0(c~ 4 ). 
If we choose Rq > such that f°(x,p) = for |x| > Rq, then 



-(cfc) W / O(o i )l B ( 0)Ro )(x + z)dpds(z) = (ct 1 B (o,r ){x + ctu) ds{u) J 0(c 4 ) 

</|z|=rf </|p|<P V y|w|=l / 

= 0(c" 4 ) 

by [23 Lemma 1], uniformly in x G M 3 , i G [0, T] and c > 1. Therefore we arrive at 

<j>x,BD{t,x) 

= -c~ 2 (cty l J J z(l-c-\z-p) + c- 2 ((z-p) 2 -l/2p 2 )+c~ 3 ((z-p)p 2 -(z-p) 3 )^) X 

x/°(x + 2,p)dpds(2) + C( C - 6 ). (5.9) 

Concerning x a , we note that f(t,x,p) = for \x\ > i? + TP =: R x . Since, by distinguishing the 
cases \x — y\ > 1 and |x — y| < 1, 

/ \z\~ 2 1 B {o,r 1 )(x + z)dz = \x - y\~ 2 l B( o iRl) (y) dy = 0(1) 

J\z\<ct J \x—y\<ct 

uniformly in a; G M 3 , £ G [0,T], and c > 1, similar computations as before show that 



<j>x,a{x,t) 



I I (z + cT 1 [p - 2(z ■ p)z] + c~ 2 [3(z ■ pfz - 3/2p 2 z - (z ■ p)p] 

J\z\<ctJ V " " " 



c 3 [4z(z ■ p)p 2 — 4(z ■ p) 3 z + (z ■ p) 2 p — p 2 p\ 



dz 

f(i(z),x + z,p)dp— + 0(c- 6 ). 
\ z \ 

In the same manner, elementary calculations using also (|2.7jl can be carried out to get 

_ 3 f f s s . dz 



(5.10) 



(j> x ,b(t,x) ■ 

<l>x,c{t,x) = -C~ 2 



-<■ I I zS((j))f(t(z),x + z,p)dp— + 0(c-^ 

\z\<ct ■' : 



(5.11) 



z\<ct 



Z® f z + cT 1 



— 2{z ■ p)z + p 



dz 



{V x (j>)f{i{z),x + z ,p)dp^ + 0(c" 4 ), 

\z\ 

4>t.BD{x, t) = -C-\ct)- 1 [ [ (l-C- 1 (z- P ) + C- 2 ((z-p) 2 -l/2p 2 ) 

J\z\=ctJ V 

f°(x + z,p) dpds(z) + 0(c- 4 ), 
(z-p) + c~ 1 (p 2 -2(z-p))^f(i(z),x + z,p)dp^ 



(j> t ,a{x,t) = C 



\z\<ct 



(5.12) 



(5.13) 



(5.14) 
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Mt,x) = -c- 2 [ [ S(4 > )f(i(z),x + z, P )d P p^ + 0(c- 2 ), (5.15) 

J\z\<ctJ \ z \ 

<f>t, c (t,x) = -c- 1 [ [ z-V x (j)f(i(z),x + z } p)dppr + 0(c- 2 ), (5.16) 

J\z\<ctJ \ Z \ 

<Ps(t,x) = -i f f f(i(z),x + z,p)dppr + 0(c- i ). (5.17) 



z <ct 



\z\ 



Next we consider the data terms treating first the gradient term. 



■>x,D 



V x .<9t(— / (j)°(x + ctu)du) J + V x ( — i <f>\x + ctu) dJ\ 
V 4?r J\u\=\ J V 47r y|w=i / 

= dt( — / V x <p°{x + ctu))duj)) + — Vx^Hx + ctuAdu 
V 47r y|w|=i / 4?r ./|«=i 

= / + // (5.18) 
Since / 2 (0, x,p) = by (|H$ . we have / / 2 (0, dp = 0. Thus we get from (fTHjl . ljVF|) and JTJJ) 

V^V) = c- 2 V a; 02(O,x) + c- 4 V :c 04(O,x) 



c 2 




M 2 ^(1 - ^)/ (0,£ + ^p)cWz + 7^ / £<9 2 /i (0, aj + z) dz. 



Using the formulas ()5.32j) - (j5.34j) below, we calculate 

J J J ' \z\~ 2 z(l - ^)f (0,x + ctu + z,p) dpdzdu 

= ~Jr j j ^)fo(0,y,p)dp J \y-x-ctu\~ 3 (y-x-ctu)dudy 



i\ u \=i 

—An f f p 2 

| z I ~ 2 z{ 1 - — ) f (0, x + z, p) dp dz, 



c 2 



1 

2^ 



\oj\ = 1 



\z\>ct 

)2 



zd t /io(0, x + ctu + z) dz dw 
2^1 J dtfMi(0,y) j \y - x - ctu\-\y - x - ctu)dudy 

~T (z~7;(c t ) 2 \ z \~ 2 z)dtVo{0,x + z)dz+-^- zd 2 fi {0,x + z)dz. 

c 4 J\ z \>ct 3 3cH J ]zl < ct 



Therefore we get 



1 4>°{x + ctu)duj) (5.19) 



I (z- ^-\z\~ 2 z)d 2 /j (0,x + z) dz + [ zd 2 /i (0, x + z) dzX 
J\z\>ct V 3 / 3c & J\ z \< ct J 

/ [ \ z \~ 2 z(l- ^)fo(0,x + z,p)dpdz + ^— [ zd 2 /i (0, x + z) dz 
J\z\>ctJ 2c 2 2c 4 J ]z] > ct 
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2c 4 
-1 

n 2 



7^ J \ z \ 2 zd t fi (0,x + z) dz + J J z(l- ^jf o (0,x + z,p)dpds(; 



'\z\=ct 

note that several terms have canceled here. A similar calculation yields 

,2 



II=-~ 

C 



3 J \A 2 j z {pth + ~ P 2*1° ) (°» g + z > P) d P dz 
+— j / zdf/j, (0,x + z)dz--^ / |z|~ 2 z<9 t 3 /i (0,2 + z)<iz 

^ c Jb|>ct J\z\>ct 



+ A / ^ 3 /i (0,a; + z)^. (5.20) 

"J c -'\z\<ct 



We need to examine the last term in (J5.20|) more closely. Employing the Vlasov-equation and 
integration by parts twice, we have 

If t f 

zdffi (0,x + z)dz = - / z(z ■ p)<9 2 /i (0, x + z) ds(z) (5.21) 



oc ->\z\<ct oc J|z|=ci 



/ [{z-p)pdtfo(0,x + z,p)dpds(z)-— 1 r[ d t U G V x <p 2 ){^,x + z) dz. 

^ J\z\=ctJ J\z\<ct 

Using the support properties and the bounds (j2.1J) - ()2.4|) as well as [23 Lemma 1] for the integrals 
over the boundary yields 

V0 D = 0(cT 2 ). (5.22) 
Therefore, (j^Hjl . (l2~7|) . together with (lo~?Ij) - (l5~T2T) and do^l gives 

V0 = C(c- 2 ) (5.23) 

and the error estimate in (J5.12)) is improved to 



(f) x , c (t,x) = -c 



-2 



! !z® 


(z + c- 1 


— 2(z ■ p)z + p 


l\z\<ct J 







dz 

{V x cf>)f{t{z),x + Z}P )dp— + 0(c- 6 ). (5.24) 

\z\ 



We next claim that also 

d t( j> = 0(c- 2 ), 
which gives the improved error estimate in (j5.11|) 

(f> x>b (t,x) = -c- 3 [ [ zS(<P)f(i(z),x + z,p)dpp- + 0(c- 6 ). (5.25) 

J\z\<ctJ \ z \ 

Combining (l5~£j) - (l5~TUj) . (|5~23J) . (|5^4l and (f53D|) - (f5^T|l . we obtain the announced formulas (|3~5|)- 
(EZD- 

In order to prove the claim and to establish representation formulas we need to examine <pr, 
and d t (pD more closely. Although we only want to give an approximation up to order c~ 3 , we have 
to take into account the terms coming from initial data of order c -4 , too, because the associated 
homogenous fields have contributions of lower order at least in the case of dt4>D- Therefore we 
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have to give the full description of 4>d an d dt'pD- Since the calculations are similar to that already 
carried out, we restrict ourself to give the formulas only: 



C J\z\>ct V ' C J\z\=ctJ 

krV/o - (z ■ p)d t f o)(0, x + z,p) dp dz 

/ / \z\~ 2 {z ■ p)d t f (0,x + z,p)dpdz 
2c 2 J\ z \ >et J 



1 

+ 2? 



--,2 



+4 / / (~dJ 2 + ?-d t f - l -^(z-p)d 2 f )(0,x + z,p)dpdz 

c J\z\<ctJ V ^ O / 



c 1 

2^ 
f3 



|z|>ct 



z >ct 



\A 1 (d t f 2 + Y d tfo^(0,x + z,p)dpdz 
\z.p)% m x + z,p)dpdz 
/ / \z\' 2 (z-p)d 2 f (0,x + z,p)dpdz. (5.26) 

OC J\z\>ctJ 



Taking into account the by now well known support properties and bounds from (|2.1j) - (|2.4jl as 
well as e.g. 



c 2 



/ \z\ 2 1b(o,Mo)( x + z ) dz < -J / 1b(o,m )(^ + z ) dz < Mc 4 > 

J\z\>ct c J\z\>ct 

^1 \z\1 B (o,m )(x + z) dz < I 1b(o,m )( x + z ) dz ^ Mc ~ 4 

c J\z\<ct c J\z\<ct 

ASHED gives 

<f> D (t,x) = — \ [ \z\~ 1 (^fi + td t fioj(0,x + z)dz + -^ I I (z-p)f Q (0,x + z,p) dpdz + G( 

(5.27) 

Combining (|5.6)L f)5 . 1 7j) and ()5.27|) we obtain 

= 0ext + 0int + 0bd + 0{c~ A ) (5.28) 



c 4 ) 



with 



<t>*tt(t,x) = — o / M 1 + #W) (0, a; + z) dz, 
0mt(*,aO = -4 / / f(i(z),x + z,p)dp d ~ 



° 2 -l'z <ct J I ; 



0bd(t,x) = -4/ [ (z ■ p)f (0,x + z,p)dpdz. 

J\z\=ct J 



'\z\=ct ■ 

Differentiating (|5.2(jjl w.r.t i and applying similar estimates yields 



d t (f) D (t,x) = 4l / ii (0,x+z)ds(z) + — [ I \z\ 2 (z •p)f (0,x + z,p)dpdz 

Ct-I\z\=ct C J\z\>ct. 
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_ i/ J( z 'P^fo(0,x + z,p)dpds(z)-^ J J p 2 j o(0,x + z,p) dp ds(z) 
+4 / / \z\- 2 (z ■ p)d t f (0,x + z,p) dpdz + 0(c~ 4 ). (5.29) 

C J\z\>ctJ 

Combining flEIJ), (jBTT^l - (fBTTB|l and (I5~2"H1) we have 

•9*0 = &,«t + 0t,iBt + ^,bd + 0(c~ 4 ) (5.30) 



with 



't.int 



(t,x) = — [ [ \z\ 2 (z ■ p){fo + td t f ){0,x + z,p)dpdz, 

C Z J\ z \>ctJ 

(t,x) = ((z-p)+c- 1 (p 2 -2(z-p))^)f(i(z),x + z,p)dp^ 

/ Z-' S ?x<j>f(t{z),X + Z,p)dpy-r, 

C J\z\<ct J \ Z \ 



<f>t,bd{t,x) = J J(z-p) 2 f (0,x + z,p)dpds(z). 



\z\=ct ■ 

Note, that because of (|5.22J) and 1)5.13]) — (|5.16j) and ()5.29|) we already have 

d t (t> = 0(c- 2 ), (5.31) 

which in turn gives 

(fH,b(t,x) = 0{c~ A ) 

<j>t,c(t,x) = -C- 1 [ [ Z-V x <j)f(i(z),x + Z,p)dppr + 0(C- 4 ), 

J\z\<ctJ \ z \ 

compare with ()5.15|1 and (|5.1fij) . 

5.2 Some explicit integrals 

We point out some formulas that have been used in the previous sections. For z G M 3 and r > 
an elementary calculation yields 

i I —l j f 47rr" 1 : r > \z\ 

\z — ruj\ duo = < , | , i . , . (5.32) 

, j ^ 4ff\z\ : r < \z\ 



Differentiation w.r.t. z gives 



M=i 



\z — ru\ 3 (z — rco)doj=l . , ^ 9 _ ' r> M . (5.33) 

' 47r|z|"^z : r < |z| v 7 



Similarly 

f ( Anr + 

/ \z — rtiA dm = < 

4ir\z\ + ^fr 2 \z\~ l : r <\z 





- ru; 


/ w|=l 





47rr + ^2 2 r 1 : r > \z\ 
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and thus by differentiation 



S7r z : r > \z\ 



z-ru\ 1 (z-ruj)duj = I 3r ' ' . (5.34) 

|u>|=i { AiTZ — ^-r z \z\ z z : r < \z\ 



Finally, for z e R 3 \ {0} also 

j\ z - v \-M~'vdv^- z (5.35) 

can be computed. 



Acknowledgments: The author is indebted to M. Kunze and H. Spohn for many discussions. 
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